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Abstract. I discuss the temperature dependence of quarkonium correlators calculated in lattice QCD. I 
show that the dominant source of the temperature dependence comes from the zero mode contribution, 
while the the temperature dependence associated with melting of bound states is quite small. I study the 
zero mode contribution quantitatively for different quark masses and show that it is well described by a 
quasi-particle model with temperature dependent heavy quark mass. As a byproduct an estimate of the 
medium dependence of the heavy quark mass is obtained. 

PACS. 12.38.Mh Quark-gluon plasma - 12.38.Gc Lattice QCD calculations 
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1 Introduction 

Quarkonium properties at finite temperature received con- 
siderable interest since the famous conjecture by Matsui 
. and Satz that color screening will lead to quarknonium 
' dissociation, which in turn can signal the onset of decon- 
, finement in heavy ion collisions [T] . This problem has been 
studied using potential models with screened tempera- 
; ture dependent potential [2,3,4,5,6,7,8,9, 10pT|[T2l[np4l 
. [T51[TBlfT71[T5] . The early studies used different phenomeno- 
' logical parametrization, while more recent ones use the 
, so-called singlet free energy of a static QQ pairs as an 
input. Since the lattice studies of the singlet free energy 
show quite strong screening (see Refs. [TOll^ l potential 
models typically predict large in-medium modification of 
quarkonium properties and/or its dissolution. 

In-medium quarkonium properties are encoded in the 
corresponding spectral functions which in turn are related 
to meson (current-current) correlation functions in Eu- 
clidean time. Euclidean correlation functions can be calcu- 
lated on the lattice. Therefore one can hope to learn some- 
thing about in-medium quarkonium properties by study- 
ing the temperature dependence of meson correlators on 
the lattice. In particular, attempts to reconstruct quarko- 
nium spectral functions using the Maximum Entropy 
Method (MEM) were presented in Refs. [2T|l22l[23l[24ll25l 
[^[77] . While at zero temperature it is possible to recon- 
struct the basic features of quarkonium spectral functions 
the application of MEM above the deconfinement temper- 
ature turns out to be problematic [26j 

The study of the spectral functions using MEM indi- 
cated that IS charmonia may survive in the deconfined 
phase till temperatures as high as (1.6 — 2.2)rc 22,23 . 
This result is closely related to the small temperature de- 
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pendence of the Euclidean correlator in the pseudo-scalar 
channel 23, 26J . On the other hand the large temperature 
dependence of the scalar and axial-vector correlators has 
been interpreted as evidence for the expected dissolution 
of the P-state charmonium. However, it has been pointed 
out that zero mode contribution, i.e. the contribution from 
the spectral functions at very small frequencies could lead 
to large temperature dependence of the Euclidean correla- 
tors [28[f29] . Therefore in order to understand quarkonium 
properties at finite temperature and their relation to col- 
lective effects, e.g. screening it is important to get a de- 
tailed understanding of the temperature dependence of the 
Euclidean correlators, including the role of the zero mode 
contribution. In this paper I present a study of the tem- 
perature dependence of quarkonium correlators calculated 
in quenched approximation on isotropic lattices. Prelimi- 
nary results of this study have been reported in conference 
proceedings t30 p 3 H [32] . The rest of the paper is organized 
as follows. I section 2 I discuss some general features of 
the Euclidean correlators and spectral functions. Section 

3 contains some details of the numerical calculations as 
well as the discussion of the temperature dependence of 
the Euclidean correlators in different channels. In section 

4 the extraction of the zero mode contribution is discussed 
and its temperature dependence is compared with a sim- 
ple quasi-particle model. Finally section 5 contains the 
conclusions. 



2 Quarkonium spectral functions and 
Euclidean correlators 

Let us consider correlation functions of meson currents 
qPg, where F — 1, 75, 7;^, 757/^, i.e. scalar, pseudo-scalar, 
vector and axial- vector currents. The relation of these cur- 
rents to different meson state is shown in Table[T] Because 
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Table 1. Quarkonium states and coefficients for the general- 
ized susceptibilities in Eq. ((2} in different channels. 
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correlation functions in Euclidean time can be considered 
as analytic continuation of the real time correlators one 
has the following relation between the Euclidean correla- 
tors and the spectral functions (e.g. see the discussion in 
Refs. PII^SPS] ) 



GHt,T) = / dua'iuj, 
Jo 



T) 



cosh(w(T - 1/2T)) 
sinh(w/2r) ■ 



(1) 



Here index i labels different channels, i = sc,ps,vc and 
ax for scalar, pseudo-scalar, vector and axial- vector re- 
spectively. In the vector and axial- vector channels we sum 
over the three spatial components. Quarkonium spectral 
functions at zero temperature contain information about 
all states with given quantum number, which contain a 
heavy quark anti-quark pair (both bound states and scat- 
tering states). At finite temperature, however, they also 
describe a scattering of the external probe (e.g. a virtual 
photon) off a heavy quark from the medium [28p33j . This 
is an analog of the Landau damping and gives a contribu- 
tion to the spectral function below the light cone {oj < k). 
In the limit of zero external momentum, — > it becomes 
X^{T)u!S{lo) in the free theory. The generalized suscepti- 
bilities X* (T) were calculated in Ref. jSSj in the free theory 
for all channels: 



x^iT) = r 

Jo 



m2 



dE„ 



Here Ep = p' 



' + m^ and npiEp) = (exp(£;p/T) + 1)" 



(2) 
The 



coefficients a' , 6* and are given in Table [T] for different 
channels. For the temporal component of the vector corre- 
lator the generalized susceptibility is just the usual quark 
number susceptibility (up to the minus sign). Interaction 
of the incoming (outgoing) quark or anti-quark with the 
medium ( multiple re-scattering) will lead to broadening 
of the delta function which becomes a Lorcntzian ^28j : 



V 



(3) 



The width of the Lorentzian ?/ ^ D^^T/m is always small 
for heavy quarks m ^ T because the diffusion constant 
D which is proportional to the mean-free path cannot be 
much smaller than the inverse temperature 1 /T |28] . The 
area under the peak at w ~ is Tx^{T) and does not 
depend on D [28j . Because the quark anti-quark pair con- 
tribute to the spectral functions only at energies uj > 2m 
quarkonium spectral functions can be written in the form 



Here (Tjj;gjj(ciJ, T) is the high energy part of the spectral 
functions which is non-zero only for uj > 2m and describes 
the propagation of bound or unbound quark anti-quark 
pairs. On the other hand al^^{uj,T) receives the domi- 
nant contribution a,t lu ~ 0. Because the width of the 
peak at w ~ is small the later gives an almost constant 
contribution to the Euclidean correlator. We can write an 
analogous decomposition for the Euclidean correlator 



G^(r,T) = GU(T,T) + GU(r,T). 



(5) 



To very good approximation GI^^{t,T) — x'(T)r, i.e. 
constant. In the next two sections I will discuss the tem- 
perature dependence of GJjjgjj and Gf^^ separately. 



3 Lattice calculations of quarkonium 
correlators 

Meson correlators have been calculated in quenched ap- 
proximation using standard Wilson gauge action and non- 
perturbatively improved Wilson action for heavy quarks. 
The analysis has been performed at three different lat- 
tices spacings and several values of the quark masses in 
the region of the charm quark mass. On the finest lattice 
also the case of bottomonium has been considered. The 
parameters used in lattice calculations, including the lat- 
tice gauge coupling (3 = and hoping parameters k are 
summarized in Table [2] The value of the clover constant 
and critical hopping parameter taken from Ref. |34j 

and interpolation. In the Table I also give the bare input 
quark mass defined as 



amo 



1 

2k 



1 



(6) 



f7X^,r) = aU(^,T) + aU(^,T). 



(4) 



At each lattice spacing meson correlators have been cal- 
culated at temperature well below the transition temper- 
ature and for several temperature values in the deconfined 
phase. The lattice volumes and the number of gauge con- 
figurations used in the analysis below Tc are also given Ta- 
ble [H The lattice spacing was set using the Sommer scale 
To = Q.5fm in contrast to Ref. "23] where it was fixed by 
the string tension. In determining the lattice spacing the 
values of rp/a calculated in Ref. [35] were used together 
with a fit to Ahon Ansatz for 6.1 < /3 < 6.92. The lattice 
spacing determined this way are given in Tabled Calcula- 
tions in the deconfined phase span the temperature range 
between l.lTc and STc. In determining the temperature 
scale T/Te I use the value roT^ = 0.7498(50) determined 
in Ref. l56] . The lattice volumes and the number of gauge 
configurations for each temperature above the transition 
temperatures are given in Table [3l 

The temperature dependence of meson correlators comes 
from the explicit temperature dependence of the spectral 
function and the trivial temperature dependence of the 
integration kernel in Eq. ([T]). Therefore to eliminate the 
trivial temperature dependence of the integration kernel 
we follow the procedure proposed in Ref. [23] and calculate 
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Table 2. Parameters used in the numerical calculations and the renormalization constants for different currents. The last 
column gives the constituent quark mass estimated in potential model (see the main text). The pseudo-scalar renormalization 
constants Zps were estimated at the scale equal to the constituent quark mass. The scalar renormalization constants Zsc were 
estimated at scales equal to half the constituent quark mass, constituent quark mass and twice the constituent quark mass 
(from left to right). 
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the so-called reconstructed correlators 



^rec ('^' 



duja {uj,T ) . , . , 7) 



where T* is some reference temperature well below Tc. Ob- 
viously if the spectral function does not change across the 
deconfinement transition the ratio /G\.^^ will be tem- 
perature independent an equal to unity. Thus we need to 
calculate the spectral functions below Tc- This has been 
done using MEM and the algorithm described in Ref. . 
The position of the first peak in the spectral function gives 
the mass of the lowest lying state in that channel. In Table 
[2] I show the mass of the ground state in the pseudo-scalar 
channel, i.e. the rjc{\S) {rn,{lS)) mass. The statistical error 
has been estimated using jackknife analysis of the peak po- 
sition. There are also systematic uncertainties in the MEM 
analysis (e.g. due to use of different priors). Therefore in 
Table [5] the combined statistical and systematic uncer- 
tainties in the determined mass values are shown For each 
input bare quark mass I estimate the constituent quark 
mass by matching the 77c(lS') {rjh{lS)) mass calculated on 
the lattice to that calculated in the potential model. A 
Cornell type potential 



V{r) = - — + crr, 
12r 



(8) 



was used with a — (1.65 — 7r/12)/rg This form gives a 
very good description of the lattice results on the static 
potential in SU(3) gauge theory. The constituent quark 
masses calculated this way are given in Table [2] 

One also has to keep in mind that meson currents 
has to be renormalized if calculated on the lattice. The 
renormalization constants for scalar, pseudo-scalar, vector 
and axial-vector currents calculated in 1-loop tadpole im- 
proved perturbation theory are given in Table [2] (see Ref. 
p3] for details of the calculations). The renormalization 
constants in the scalar and pseudo-scalar channels depend 
on the renormalization scale which was chosen to be the 
constituent heavy quark mass. In the scalar channel we 
will also need the renormalization constant estimated at 
scale equal to half and twice the constituent quark mass. 



To demonstrate the general features of the spectral 
functions below Tc in Fig. [1] I show the spectral functions 
calculated for (3 = 7.192. In addition to ground state peak 
they show a second peak which may be due to combination 
of several excited states and two broad bumps at even 
higher energy. The later are most likely artifacts of limited 
statistics. In Ref. [26], where the spectral functions have 
been calculated with significantly more statistics the high 
energy part of the spectral functions appears to be much 
smoother. 

As mentioned above the ratio G{T,T)/Grec{T,T) can 
be used to study the temperature dependence of the quarko- 
nium correlatorsEI. However, the temperature dependence 
of the correlators can be induced by the high energy part 
of the spectral functions and the low energy part of the 
spectral functions, which gives the zero mode contribu- 
tion. The zero mode contribution is absent in the deriva- 
tive of the correlator with respect to r. Therefore the best 
way to study the temperature dependence of the correla- 
tors induced by change of bound state properties and/or 
its dissolution is to consider the ratio of the derivatives of 
the correlators G"(t, r)/G'^g^(r, T). In the next two sub- 
section we are going to discuss this ratio for pseudo-scalar 
channel as well as for scalar and axial- vector channels re- 
spectively. 



3.1 The temperature dependence of the pseudo-scalar 
correlators 

The temperature dependence of the charmonium and bot- 
tomonium correlators in pseudo-scalar channel has been 
studied in Refs. [23ll26p 30]. In calculations on isotropic 
lattices the charmonium correlators do not show any tem- 
perature dependence up to temperatures 1.5Tc. More pre- 
cisely the ratio G{T,T)/Grec{T,T) is one within statisti- 
cal errors up to this temperature. At higher temperatures 
we see deviation from unity in this ratio. In particular 
at 3Tc this ratio deviates from one by about 10% f23^. In 
anisotropic lattice calculations the ratio G(t, T) / Grec{T, T) 



This value of a follows from the definition of the Sommer 



scale r'^Vr{r)\ 



=""0 



= 1.65. 



In the remainder of this section the subscript i labeling the 
channel will be omitted. It should be clear from the context 
which channel is discussed. 
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Fig. 1. The spectral functions in different channels calculated 
using 48^ X 64 lattice ai (3 = 7.192. The vector and axial- 
vector spectral functions are divided by two. The horizontal 
line shows the spectral function in the massless free theory 
a{w)/uj^ = 3/{8tv^). 



is also close to one up to temperatures 1.5Tc although sta- 
tistically significant deviation from unity are seen at the 
level of few percent. At higher temperatures, however, one 
sees significant deviations from unity which are also larger 
than those observed in isotropic lattice calculations and 
are present at smaller values of r [26]. These deviations 
are larger on coarser lattices. In the bottomonium case no 
temperature dependence was found up to STc- In Figure 
[2] I show the ratio of the derivatives G'{t,T)/G'^^^{t,T) 
calculated at /3 = 7.192. The results for this ratio cal- 
culated on anisotropic lattices with larger spatial lattice 
spacings are also shown [26| . As one can see from the fig- 
ure on isotropic lattices the ratio of the derivatives docs 
not show any temperature dependence and agrees with 
one within statistical errors. The temperature dependence 
of the pseudo-scalar correlator seen in anisotropic lattice 
calculations is also greatly reduced and the agreement be- 
tween isotropic and anisotropic calculations is is better 
for this ratio. One possible explanation of this behavior is 
the presence of a small negative zero mode in the pseudo- 
scalar channel. In the continuum there is no zero mode 
contribution in the pseudo-scalar at leading order (c.f. Ta- 
ble [1]) and this is true also on the lattice [37] . While we 
do not expect to have a negative zero mode contribution 
even at higher orders in the continuum theory, the situa- 
tion may be different on the lattice. There could be a neg- 
ative zero mode contribution in the pseudo-scalar chan- 
nel which vanishes in the continuum limit The fact that 
deviations from unity in the ratio G/Grec increases with 
increasing the lattice spacing makes this scenario a plausi- 
ble one. Apparently modifications of bound state proper- 
ties is not reflected in the temperature dependence of the 
correlators, i.e. Ghigh(''', 7") is temperature independent at 
the level of few percent. The temperature dependence of 
the pseudo-scalar correlators is consistent with potential 
model calculations [T5]ll8j. 



3.2 Scalar and axial-vector correlators 

The temperature dependence of scalar and axial-vector 
correlators has been studied in terms of G/Grec using 
isotropic [23] as well as anisotropic lattices [26]. The ra- 
tio G/Grec is temperature independent in the confined 
phase |26j and shows large enhancement in the deconfined 
phase p3ll26| . This large enhancement is present both 
in charmonium and bottomonium correlators. To elim- 
inate the zero mode contribution I have calculated the 
derivative of the correlators and the corresponding ratio 
G'{T,T)/G'rec{'T,T). The numerical results for this ratio 
at /? = 7.192 are shown in Figure |3l Here I also show 
the results from anisotropic lattices [53]. There is good 
agreement between the results obtained from isotropic 
and anisotropic lattices. As one can see from the Figure 
G'{T,T)/G're^{T,T) shows very little temperature depen- 
dence and is close to unity. This means that almost the 
entire temperature dependence of the scalar and axial- 
vector correlators is due to zero mode contribution and 
GhighlT, T) is temperature independent. One may won- 
der if the temperature independence of Ghigh (t, T) implies 
the survival of the P-state quarkonium up to temperature 
as high as 3Tc. The P-state quarkonium has binding en- 
ergy that is significantly smaller than that of the S-state. 
Therefore we expect it to melt at lower temperatures. In- 
deed, potential model calculations show that the disso- 
ciation temperatures for the P-states is significant lower 
than for S-states and none of the potential models can 
get dissociation temperatures as high as 3Tc for P-states. 
The potential model calculations [2l[3l[4l[5l[6l[7K8l[9l[T0pT] 
12,13,14,15,16,17] take into account the temperature de- 
pendence of the real part of the potential. At finite temper- 
ature the potential has also an imaginary part [38ll39] . The 
imaginary part of the potential further weakens the bound 
state contribution '40' . The very weak temperature depen- 
dence of the quarkonium correlators despite of melting of 
the bound states looks puzzling at a first glance. However, 
the calculation of quarkonium correlators within potential 
model shows that the high energy part of the correlator 
almost does not change with the temperature even when 
the P-state quarkonium melts \\^ . 

Note that the temperature independence of meson cor- 
relators is true only for heavy quarks. In the light quark 
sector meson correlators show significant temperature de- 
pendence |42P43). 

4 Estimating the zero mode contribution 

As discussed in the previous section the high energy part 
of the quarkonium correlators turns out to be temperature 
independent to very good approximation. Therefore we 
could use the spectral functions calculated below Tc to es- 
timate it, i.e. we can assume that G\^:^^^{t, T) ~ G^^^It, T). 
Then the low energy part of the correlators, i.e. the zero 
mode contribution can be evaluated as G\^^{t^ T) = G'^{t, T) — 
G\.^^{t,T). It is important to estimate the systematic er- 
rors in this difference due to inaccuracies in the deter- 
mination of the spectral functions in the confined phase. 
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Fig. 2. The ratio of the derivatives G'(t, T)/Grec(T, T) in the 
pseudo-scalar channel calculated on isotropic lattices at /3 = 
7.192 (top) as well results from anisotropic lattice calculations 
at /3 = 6.5 [26] (bottom). 
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Fig. 3. The ratio of the derivatives G' {t,T) /G'r^c{r,T) in the 
scalar channel (top) and axial-vector channel (bottom) calcu- 
lated at /3 = 7.192. The results from anisotropic lattice calcu- 
lations at /3 = 6.5 [26] are also shown (open symbols). 



Table 3. Finite temperature lattices and the corresponding 
number of gauge configurations used in the analysis of quarko- 
nium correlators in the deconfined phase. For each lattice vol- 
ume the corresponding value of the temperature is also shown. 





lattice 


#config. 




6.499 


48^ X 12 


50 


1.23 




48^ X 10 


46 


1.48 
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80 


1.47 
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100 
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48^ X 12 


90 
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i.e. GI.^^{t,T), as well as due to the small temperature 
dependence of (Thigh (w). To estimate the errors in the re- 
constructed correlator I compared several spectral func- 
tions calculated using different Bayesian priors and differ- 
ent jackknife sub-samples. The uncertainties in G\.^^{t, T) 
turn out to be smaller than 1%. To estimate the uncer- 
tainties due to temperature dependence of (Thigh(w), i.e. 
due to the fact that it is slightly different from the zero 



temperature spectral functions we consider the pseudo- 
scalar channel. Here we do not expect to have a zero 
mode contribution in the continuum limit and therefore 
G*(t,T) - Gi^^{T,T) ~ 0. At 1.5Tc this difference is in- 
deed zero within statistical errors and we used it to es- 
timate the systematic errors in other channels. At higher 
temperatures as discussed in the previous section we may 
have a negative zero mode contribution in the pseudo- 
scalar channel at finite lattice spacing. But this zero mode 
contribution is much smaller than the zero mode contri- 
bution in other channels. 

In Figure [3] I show the estimated zero mode contribu- 
tion GI^^{t,T)/T^ for the vector and axial- vector chan- 
nels. In the figure both statistical (thin error bars) and sys- 
tematic (thick error bars) uncertainties are shown. Both 
uncertainties increase with decreasing r. This is due to 
the fact that the relative contribution of the zero mode 
decreases fast at small Euclidean times. 

To the first approximation Gl^^ (t, T) /T^ should be 
constant. We see that the numerical results for this quan- 
tity are indeed compatible with constant behavior in t 
within the estimated errors. The curvature of G\^^{t,T) 



6 



Peter Petreczky: On temperature dependence of quarkonium correlators 



0.05 
0.045 - 

0.04 - 
0.035 

0.03 - 
0.025 



D=1/(2 7iT) 




0.25 0.3 0.35 0.4 0.45 0.5 



0.18 
0.17 
0.16 
0.15 
0.14 
0.13 
0.12 
0.11 
0.1 



Gfow (^)/t' 



0.25 0.3 0.35 0.4 0.45 



0.5 



Fig. 4. The low energy part of the vector (top) and axial- 
vector (bottom) correlators as function of r. The thin error 
bars correspond to statistical error, while the thick error bars 
to systematic errors (see the main text). 



is controlled by the diffusion constant D 28J. However, 
even for the smallest possible value of the diffusion con- 
stant D — 1/{2ttT) (one thermal wavelength) the cur- 
vature of G\^^{t,T) is quite small as can be seen from 
the Figure Si Current lattice data are clearly not precise 
enough to constrain the value of D, they can only tell 
about the area under the peak at w ~ 0. Therefore we 
could take the value of low energy part of the correlator 
at the midpoint as an estimate for the zero mode contri- 
bution, GI^{t = l/(2r),T) = Tx'(T). In the following 
I will discuss the temperature dependence of x*(r) for 
different quark masses and in different channels. 

Since the zero mode contribution is related to the prop- 
agation of single (unbound) heavy quark in the medium 
it is natural to ask to what extent its temperature depen- 
dence can be described in terms of a quasi-particle model. 
For this purpose let us consider the temporal component of 
the vector correlator G"'^°. In this case there is no high en- 
ergy component and 0""° = -Tx^^^iT) = -Tx{T), with 
x(T) being the heavy quark number susceptibility. In the 
free theory x/T'^ depends only on m/T. Thus matching 
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Fig. 5. The quark number susceptibility as function of 
T /rricsiT) (top) and the thermal mass correction (bottom) as 
function of the temperature. The solid black line corresponds 
to the quasi-particle model. The dashed line and the band cor- 
respond to the perturbative prediction of the thermal mass 
correction. 



the lattice data on C"^", we can define and effective tem- 
perature dependent heavy quark mass mcfi(T). It is also 
convenient to define an effective thermal mass correction 
(5meff(r) = mcff(T) — m. In Figure [5] I show the quark 
number susceptibility as function of T/rUee as well as the 
effective thermal mass correction Sm(.f{{T). The thermal 
correction to the heavy quark mass is largest at Tc and 
monotonically decreases with increasing temperature. It 
also depends on the heavy quark mass, namely it decreases 
with increasing quark mass. For very heavy quarks the 
thermal mass correction is given by [38 l l39 1 l44] 



c 4g2(T) 

dmcff = -7;—, mo, 

3 47r 



(9) 



the free theory expression for x 



vcO 



X{T) (Eq. m) to 



with Tod = g{T)T being the perturbative Debye mass 
{Nf = because we work in quenched approximation) . It 
is interesting to compare this leading order perturbative 
prediction with the numerical results on 5mcti{T). For this 
we have to determine the running coupling g{T). I used the 
2-loop running coupling evaluated at scale /i = 27rT as well 
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Fig. 6. The thermal velocity of heavy quarks (top) and 
the zero mode contribution to the ajdal-vector correlator 
Gfow/r^ = X^^/T^ (bottom). The lines show the prediction 
of the quasi-particle model with mcff(T). The open symbols 
show the thermal velocity squared estimated on anisotropic 
lattices 1261. 



as the relation Tc/Ajjg = 1.15(5) [IS]. The corresponding 
perturbative predictions are shown as the dashed Une in 
Figure O To estimate the uncertainties in this perturba- 
tive prediction I also evaluated the coupling at /i = ttT 
and fx = AttT . This is shown as the band in Figure [5l 

Now I discuss the temperature dependence of the zero 
mode contribution in the vector and axial- vector channels. 
In general the zero mode contribution depends on both the 
temperature and the quark mass. In quasi-particle picture, 
however, it is a function of T/rricft only. Therefore in Fig- 
ure[6]I show the numerical results for x"^ (T) as function of 
T/rricS- Here I used the temperature dependent effective 
quark masses determined above. We see that x''^(r)/T^ 
is indeed function of T/nieff only and agrees very well 
with the prediction of the quasi-particle model shown in 
the Figure as the black line. 

Similar analysis has been done in the vector channel. 
In this case I analyzed the ratio G\^^ (T) / G'"'^^ which has 
a simple physical interpretation as the averaged thermal 
velocity squared of the heavy quark u^j^ (T) . This is because 



due to the heavy quark mass the Fermi-Dirac distribution 
in Eq. ([2]) can be replaced by the Boltzmann distribution 
and we have 



p 



I \ d'pe 



-Ep/T 



(10) 

The numerical results for v^,^{T) = G'(;;^{T) / G"""^ {T) are 
shown in Figure [6l As one can see also v^^iT) is a function 
of T/rrics only and agrees with the prediction of a quasi- 
particle model shown as the black curve. The zero mode 
contribution in the vector channel has also been studied 
on anisotropic lattices [15] and in Figure [H] I also show 
the corresponding results for thermal velocity squared. To 
estimate rrioff in this case I used 1.19GeV for the con- 
stituent charm quark mass needed to get the charmonium 
spectrum right in the potential model [15] and SniesiT) 
determined above. For T ^ m we have v'^^^ T/m but as 
one can see from the figure it is not a good approximation 
even for the b-quark. 

Finally I discuss the zero mode contribution in the 
scalar channel. Unlike the vector and axial-vector corre- 
lators the scalar correlator depends on the choice of the 
renormalization scale. It is not obvious what is the suit- 
able choice of the scale when comparing the zero mode 
contribution to the prediction of the quasi-particle model. 
Since the zero mode contribution is related to the propa- 
gation of a single quark the choice fi = m for the renor- 
malization scale seems to be a natural one. To get an idea 
about the sensitivity of the conclusions to the choice of the 
renormalization scale I have also considered fi = m/2 and 
fi = 2m. Using the renormalization constants given in Ta- 
ble[2] the zero mode contribution in the scalar channel has 
been evaluated at scales fi = m/2, to, 2m. The numerical 
results for the zero mode contribution to the scalar cor- 
relator as function of T/mcs are shown in Figure [Tj We 
see again that the zero mode contribution is function of 
T/TOcff only and there is a reasonable agreement with the 
prediction of the quasi-particle model shown as a black 
line. 

One way to quantify the zero mode contribution to the 
scalar correlator which is independent of the choice of the 
renormalization scale is to consider the ratio Gf^^/G^g^ 
evaluated at r = 1/(2T). The numerical results for this 
ratio are shown in Figure [H As one can see from the figure 
this quantity decreases with increasing temperature and 
decreasing quark masses. There is no large dependence 
lattice spacing in this quantity. In the Figure I also show 
the results obtained on anisotropic lattices at /3 = 6.5 j26j . 
which are significantly larger than the results obtained on 
isotropic lattice. This is most likely due to the fact that 
the spatial lattice spacing at (3 = 6.5 is too coarse for 
studying the zero mode contribution. 



5 Conclusions 

In this paper the temperature dependence of quarkonium 
correlation function has been discussed. It has been shown 
that the temperature dependence of the high energy part 
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Fig. 7. The zero mode contribution to the scalar correlator in 
units of T'^ as function of T/nicsiT). The black line is the pre- 
diction of quasi-particle model. The filled symbols correspond 
to the scalar correlator evaluated at scale /j, = m, while the 
open symbols correspond to the scalar correlator evaluated at 
scale n = m/2 and n = 2m. 
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Fig. 8. The ratio of the zero mode contribution of the scalar 
correlator to the corresponding reconstructed correlator at r = 
l/(2r). Open symbols correspond to the results obtained on 
anisotropic lattices [26j . 



of the spectral function, for example, the melting of reso- 
nances does not lead to a large change in the correlation 
function. The dominant source of the temperature depen- 
dence of quarkonium correlators is the zero mode contri- 
bution. This contribution has been studied quantitatively 
on the lattice for the first time. In general, it is expected 
that this contribution depends on the temperature and 
the quark mass. We have found, however, that it is the 
function of rrics /T only and is well described by a quasi- 
particle model down to temperatures as low as l.lTc. Since 
the the quasi-particle mode is so successful in describing 
the zero mode contribution to the quarkonium correlators 
it would be of great interest to calculate it systematically 
in improved perturbation theory [46j . 
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